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Abstract 

In future collider experiments at the TeV scale, large logarithmic corrections 
originating from massive boson exchange can lead to significant corrections to 
observable cross sections. Recently double logarithms of the Sudakov-type were 
resummed for spontaneously broken gauge theories and found to exponentiate. 
In this paper we use the virtual contributions to the Altarelli-Parisi splitting 
functions to obtain the next to leading order kernel of the infrared evolution 
equation in the fixed angle scattering regime at high energies where particle 
masses can be neglected. In this regime the virtual corrections can be described 
by a generalized renormalization group equation with infrared singular anoma- 
lous dimensions. The results are valid for virtual electroweak corrections to 
fermions and transversely polarized vector bosons with an arbitrary number of 
external lines. The subleading terms are found to exponentiate as well and are 
related to external lines, allowing for a probabilistic interpretation in the mass- 
less limit. For Z-boson and 7 final states our approach leads to exponentiation 
with respect to each amplitude containing the fields of the unbroken theory. For 
longitudinal degrees of freedom it is shown that the equivalence theorem can be 
used to obtain the correct double logarithmic asymptotics. At the subleading 
level, corrections to the would be Goldstone bosons contribute which should be 
considered separately. Explicit comparisons with existing one loop calculations 
are made. 
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1 Introduction 



Future high energy experiments at the LHC and possibly a hnear colhder (TESLA 
or the NLC for instance), will probe the full non-Abelian nature of the electroweak 
Standard Model (SM). Thus one has to view the photon in particular as a particle 
with non-Abelian character. At energies much larger than the weak scale, Sudakov 
logarithms, originating from vector boson exchange, can lead to significant radiative 
corrections. The double logarithms (DL) can be of order (9(20 %) at one loop in the 
TeV range and a few % at the two loop level. In addition, subleading corrections can 
also be significant, especially if the experimental accuracy is of the order of 0{1 %). 

As of this writing, there is no complete two loop calculation in the electroweak 
theory due to the complexity of the number and nature of processes involved. It is 
thus of considerable interest to investigate terms which are potentially large and which 
can be resummed to all orders. In Ref. the leading DL corrections were calculated 
and found to exponentiate. The results were obtained by using the infrared evolution 
equation method calculated with the massless fields of the unbroken theory. The 
equation has a different kernel depending on the value of the infrared cutoff. 

There are, however, some important differences of the electroweak theory with 
respect to an unbroken gauge theory. Since the physical cutoff of the massive gauge 
bosons is the weak scale M = M^, pure virtual correction lead to physical cross 
sections depending on the infrared "cutoff". Only the photon needs to be treated 
in a semi-inclusive way. Additional complications arise due to the mixing involved 
to make the mass eigenstates and the fact that at high energies, the longitudinal 
degrees of freedom are not suppressed. Furthermore, since the asymptotic states are 
not group singlets, it is expected that fully inclusive cross sections contain Bloch- 
Nordsieck violating electroweak corrections 

In this paper we extend the method of Ref. [|^ to the next to leading order for 
the case of virtual corrections at high energies where we can neglect particle masses. 
In addition, we show that the results of Ref. are also valid for longitudinal de- 
grees of freedom, which at first sight, is far from obvious. The connection between the 
calculation performed in the massless theory and the longitudinal degrees of freedom 
is provided by the Goldstone boson equivalence theorem. Another complication in 
comparison with the unbroken non-Abelian case is the mixing of the mass-eigenstates. 
Especially for the external Z-boson and the photon states, the corresponding correc- 
tions in general don't factorize with respect to the original amplitude. We indicate 
how the corrections are given to subleading level in terms of the fields of the unbroken 
theory. 

Finally we compare our results with existing one-loop corrections in the high en- 
ergy approximation for on-shell VT-pair production in electron-positron scattering. Al- 
though this comparison constitutes a strong test of our approach it should be mentioned 
that it would be extremely helpful to compare the results obtained in this approach 
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with a general method in terms of the physical fields. In this context a two loop DL- 
calculation would help clarify the situation with contradicting results in the literature 
and a subleading one loop approach would give further support to the results 
presented in this work. 

2 Logarithmic corrections in non-Abelian theories 

In this section we are concerned with virtual double and single logarithmic corrections 
to scattering amplitudes in massless non-Abelian theories at fixed angle with all invari- 
ants large with respect to an infrared cutoff /i, i.e. fi"^ <C Sj^i = 2pjPi ~ s. It must be 
emphasized that in high energy collider experiments there are also contributions de- 
pending on angular variables (i.e. u/t etc.) which can be of genuine subleading nature 
0. The philosophy adopted here is that terms of the type log ^ log j etc. should be 
calculated exactly at least at the one loop level. For the higher order terms below, we 
are only concerned with the log^ behavior with /i^ <^ s ~ |t| ~ \u\. All mass terms 
are neglected, i.e. we assume < fi. 

We begin by reviewing the general method for virtual corrections in the DL- 
approximation following the approach of Ref . . 

2.1 Double logarithmic corrections 

Sudakov effects have been widely discussed for non-Abelian gauge theories, such as 
SU{N) and can be calculated in various ways (see, for instance, 0). A general method 
of finding the DL asymptotics (not only of the Sudakov type) is based on the infrared 
evolution equations describing the dependence of the amplitudes on the infrared cutoff 
/i of the virtual particle transverse momenta [@]. This cutoff plays the same role as 
the fictitious photon mass A in QED, but, unlike A, it is not necessary that it vanishes 
and it may take an arbitrary value. It can be introduced in a gauge invariant way by 
working, for instance, in a finite phase space volume in the transverse direction with 
linear size / ~ Instead of calculating asymptotics of particular Feynman diagrams 
and summing these asymptotics for a process with n external lines it is convenient 
to extract the virtual particle with the smallest value of in such a way, that the 
transverse momenta of the other virtual particles are much bigger 

k'l > fci > (1) 

For the other particles plays the role of the initial infrared cut-off fi'^. 

In particular, the Sudakov DL corrections are related to the exchange of soft gauge 
bosons. For this case the integral over the momentum k of the soft (i.e. \k^\ -C 
virtual boson with the smallest k± can be factored off, which leads to the following 
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infrared evolution equation: 



2 (27r)'^ j i^.j^iJs-»k\-»p,^ + «e {kpj){kpi) 
xT"(j)T'^(0-M(pi,...,P„;fcl), (2) 

where the amphtude ...,Pn] ^1) on the right hand side is to be taken on the mass 

shell, but with the substituted infrared cutoff: fi^ — > k^. The generator T"'{l){a = 
1, A^) acts on the color indices of the particle with momentum pi. The non-Abelian 
gauge coupling is g. In Eq. (||), and below, k± denotes the component of the gauge 
boson momentum k transverse to the particle emitting this boson. It can be expressed 
in invariant form as fc^ = mm{(kpi){kpj) / (pipj)) for all j ^ I. 

The above factorization is related to a non-Abelian generalization of the Gribov 
theorem^ for the amplitude of the Bremsstrahlung of a photon with small transverse 
momentum k± in high energy hadron scattering [Q. 

The form in which we present Eq. (0) corresponds to a covariant gauge for the gluon 
with momentum k. Formally this expression can be written in a gauge invariant way 
if we include in the sum the term with j = I (which does not give a DL contribution). 
Indeed, in this case we can substitute piPj by —pfp'^dfj,u{k), where the polarization 
matrices of the boson dfj^^ik) in the various gauges differ by the terms proportional 
to k^ or k^ giving a vanishing contribution due to the conservation of the total color 
charge J2aT°' = 0. Thus we have the possibility of choosing appropriate gauges for 
each kinematical region of quasi-coUinearity of k and pi. We can, however, use (0) as 
well, noting that in this region for j ^ I we have PjPi/kpj ~ Ei/u, where Ei is the 
energy of the particle with momentum pi and a; the frequency of the emitted gauge 
boson, so that: 

M{pu...,Pn;^^') = MBorniPu -^Pn) - -r^Y. T^/ " 

X CiM{pi,...,Pn;kl) , (3) 

where Ci is the eigenvalue of the Casimir operator T°-{l)T"-{l) {Ci = Ca for gauge 
bosons in the adjoint representation of the gauge group SU{N) and Ci = Cp for 
fermions in the fundamental representation). In this last step we also used the identity 
YJj=iT"'{j)M.{pii ■■■■iPn] k\) = 0, corresponding to the conservation of the total group 
charge. The integral over d'^k was written in terms of the Sudakov components accord- 



ing to the discussion in section |2.2| upon replacing the longitudinal component u with 
the boson on-shell expression suv = /i^ + fc^. Thus, in Sudakov DL corrections there 
are no interference effects, so that we can talk about the emission (and absorption) of 

^The non-Abelian generalization of Gribov's theorem is given in Ref. [0, together with a descrip- 
tion of its essential content. 
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a gauge boson by a definite (external) particle, namely by a particle with momentum 
almost coUinear to k. 

The differential form of the infrared evolution equation follows immediately from 

i): 

0^—^^ = K{^i)M{pi,...,Pn]^). (4) 

where 
with 

^K^'/^') = (|)^^'l°g';^- (6) 

PVi is the probability to emit a soft and almost coUinear gauge boson from the particle 
/, subject to the infrared cut-off /i on the transverse momentum Note again that 
the cut-off /i is not taken to zero. To logarithmic accuracy, we obtain directly from 

i): 

The infrared evolution equation (^) should be solved with an appropriate initial con- 
dition. In the case of large scattering angles, if we choose the cut-off to be the large 
scale s then clearly there are no Sudakov corrections. The initial condition is therefore 

M{pi, ...,Pn; S) = MBorniPl, ■■■,Pn), (8) 

and the solution of (Q) is thus given by the product of the Born amplitude and the 
Sudakov form factors: 

= A^Born(Pl,...,Pn)exp (9) 

Therefore we obtain an exactly analogous Sudakov exponentiation for the gauge group 
SU{N) to that for the Abehan case g. 



2.2 Soft divergences in the massless theory 

In this section we briefly review the types of soft, i.e. -C i/i, divergences in loop 
corrections with massless particles. In general, those contributions, unlike the coUinear 
logarithms, can be obtained by setting all k dependent terms in the numerator of tensor 
integrals to zero (since the terms left are of the order of the hard scale s). Thus it is 
clear that the tensor structure which emerges is that of the inner scattering amplitude 
in Fig. D taken on the mass-shell, times a scalar function of the given loop correction. 
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In the Feynman gauge, for instance, we find for the well known vertex corrections 
the familiar three-point function Cq and for higher point functions we note that in the 
considered case all infrared divergent scalar integrals reduce to Cq multiplied by factors 
of ^ etc.. The only infrared divergent three point function is given by 

Co(s//^ ) = y^^^^ (27r)4 + ie^k'' + 2pjk + ie){k'' - 2pik + ie) 

It is now convenient to use the Sudakov parametrization for the exchanged virtual 
boson: 

k = vpj +upi + k± (11) 
For the boson propagator we use the identity 

7rd{suv-kl) (12) 



SUV — kj_ + ie SUV — kj^ 

writing it in form of the real and imaginary parts (the principle value is indicated by 
V). The latter does not contribute to the DL asymptotics and at higher orders gives 
subsubleading contributions. Rewriting the measure as d'^k — (Pk±(Pk\\ with 

d^k± = \k±\d\k±\d(i) = l-dkld(p = TTdkl (13) 
d\\ = \d{k^,k'')ld{u,v)\dudv^^-dudv (14) 

where we turn the coordinate system such that the Pj, plane corresponds to 0, x and 
the z coordinates to the k^ direction so that it is purely spacelike. The last equation 
follows from = 0, i.e. pf^ f=:i pf^, and 

(PjoP/. - Vi^Vj.f ~ (PioP/o - ViJPj.f = (PjPiT = «/2 (15) 

The function Cq{s/ jj?) is fastly converging for large k^ and we are interested here in 
the region /i^ <^ s in order to obtain large logarithms. Then logarithmic corrections 
come from the region fc^ <^ s\v\ -C s (the strong inequalities give DL, the simple 
inequalities single ones) and we can write to logarithmic accuracy: 



roo roo poo 

Co(s//^') = 7^77^ / du dv dklx 

Z[Z7r) J-oo J-oo J 



{suv — + ie){suv — k^-\- su -\- ie){suv — kj_ — sv + ie) 
si^ r^dur^d^ ^ _ _ 

inr J-oo SU J-oo SV J-oo 



2(27r)^ 

i du dv^, 

— —Q{suv - 11^) 
xrs J-i u J-i V 



2(47r) 
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1 du dv 



(AttYs Jli^ U Jii^ V 

^ / s su 



I 



W 4 (16) 



2(47r)2s ^ /i' 

Thus, no single soft logarithmic corrections are present in Cq{s/h'^). In order to see 
that this result is not just a consequence of our regulator, we repeat the calculation 
for a fictitious gluon mas^. In this case we have 

Co{s/X ) = J (27r)4 (P - A2 + ie)(P + 2pjk + ie){k^ - 2pik + ie) ^ ' 

It is clear that Co(s/A^) contains soft and coUinear divergences [k \\ pj^i) and is regu- 
lated with the cutoff A, which plays the role of fi in this case. Integrating over Feynman 
parameters we find: 

We are only interested here in the real part of loop corrections of scattering amplitudes 
since they are multiplied by the Born amplitude and the imaginary pieces contribute to 
cross sections at the next to next to leading level as mentioned above. In fact, the minus 
sign inside the double logarithm corresponds precisely to the omitted principle value 
contribution of Eq. (0) in the previous calculation. Thus, no single soft logarithmic 
correction is present in the case when particle masses can be neglected. 

This feature prevails to higher orders as well since it has been shown that also in 
non-Abelian gauge theories the one- loop Sudakov form factor exponentiates 0. 

In case we would keep mass-terms, even two point functions, which in our scheme 
can only yield coUinear logarithms, would contain a soft logarithm due to the mass- 
renormalization which introduces a derivative contribution |T^. In conclusion, all 
leading soft corrections are contained in double logarithms (soft and coUinear) and 
subleading logarithmic corrections in a massless theory, with all invariants large {sj^i = 
'^PjPi ~ ^{^)) compared to the infrared cutoff, are of the coUinear type or renormal- 
ization group logarithms. 



^Note that this regulator spoils gauge invariance and leads to possible inconsistencies at higher 
orders. Great care must be taken for instance when a three gluon vertex is regulated inside a loop 
integral. 
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Figure 1: In an axial gauge, all collinear logarithms come from corrections to a partic- 
ular external line (depending on the choice of the four vector satisfying n'^Al = 0) 
as illustrated in the figure. In a covariant gauge, the sum over all possible insertions 
is reduced to a sum over all n-external legs due to Ward identities. Overall, these 
corrections factorize with respect to the Born amplitude. 

2.3 Virtual logarithmic corrections from the Altarelli-Parisi 
splitting functions 

In an axial gauge, collinear logarithms are related to corrections on a particular ex- 
ternal leg depending on the choice of the four vector riu [|1T|. A typical diagram is 



depicted in Fig. |I[ In a general covariant gauge this corresponds (using Ward iden- 
tities) to a sum over insertions in all n external legs We can therefore adopt the 
strategy to extract the gauge invariant contribution from the external line corrections 
on the invariant matrix element at the subleading level. The results of the previous 
section are thus important in that they allow the use of the Altarelli-Parisi approach 
to calculate the subleading contribution to the evolution kernel of Eq. (H). We are 
here only concerned with virtual corrections and use the universality of the splitting 
functions to calculate the subleading terms. For this purpose we use the virtual quark 
and gluon contributions to the splitting functions P^q{z) and P^g{z) describing the 
probability to emit a soft and/or collinear virtual particle with energy fraction z of the 
original external line four momentum. The infinite momentum frame corresponds to 
the Sudakov parametrization with lightlike vectors. In general, the splitting functions 
Pba describe the probability of finding a particle B inside a particle A with fraction z 
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of the longitudinal momentum of A with probability Vba to first order 



dVBA{z) = -^PsAdt (19) 

ZTT 



where the variable t = log ^ for our purposes. It then follows that 



dVBA{z) = — 2^ p dlogk^ (20) 

spins -L 



where Va ,b+c denotes the elementary vertices and 



|2 



p,,u) = £(i_i)^n<i=£±^ (21) 

spins 1. 

The upper bound on the integral over dk^ in Eq. ( pil| ) is s and it is thus directly related 
to dt. Regulating the virtual infrared divergences with the transverse momentum cutoff 
as described above, we find the virtual contributions to the splitting functions for 
external quark and gluon lines: 

P^^iz) = (^-2 log -^ + 3^ 5(1-^) (22) 
= C^(-21og-^ + ^/?o'^^°)5(l-z) (23) 

The functions can be calculated directly from loop corrections to the elementary pro- 
|T3], |T^, |T^ and the logarithmic term corresponds to the leading kernel of section 



cesses 



2.1. We introduce virtual distribution functions which include only the effects of loop 



computations. These fulfill the Altarelli-Parisi equations^ 

MPi'M (24) 

/y,t)P^siy) (25) 

The splitting functions are related by Pba = Pba + ^ba^ where R denotes the con- 
tribution from real gauge boson emissionR. Pba is free of logarithmic corrections and 



dq{z,t) 


_ as 


f^dy 


dt 


2tt 


Jz y 


dg{z,t) 


as 


dy 


dt 


27r 


Jz y 



■^Note that the off diagonal sphtting functions Pqg and Pgg do not contribute to the virtual proba- 
bilities to the order we are working here. In fact, for virtual corrections there is no need to introduce 
off-diagonal terms as the corrections factorize with respect to the Born amplitude. The normalization 
of the Eqs. ( |2^ ) and ( p3| ) corresponds to calculations in two to two processes on the cross section 
level with the gluon symmetry factor i included. The results, properly normalized, are process 
independent. 

^P„„ was first calculated by V.N. Gribov and L.N. Lipatov in the context of QED nM. 
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Figure 2: Feynman diagrams contributing to the infrared evolution equation ( ^2]) for 
a process with n external legs. In a general covariant gauge the virtual gluon with 
the smallest value of fc^ is attached to different external lines. The inner scattering 
amplitude is assumed to be on the mass shell. 



positive definite. The subleading term in Eq. ( p3D indicates that the only subleading 
corrections in the pure glue sector are related to a shift in the scale of the coupling. 
These corrections enter with a different sign compared to the conventional running 
coupling effects. The renormalizations with respect to the Born amplitude as well as 
the ones belonging to the next to leading terms at higher orders will be indicated be- 



low by writing Q;s(/i^ 



a. 



[s). For fermion lines there is an additional subleading 



correction from collinear terms which is not related to a change in the scale of the 
coupling. 

Inserting the virtual probabilities of Eqs. (22) and (|23|) into the Eqs. (24) and (25) 
we find: 



g(i,t) 



goexp 
5(0 exp 



log ^ - 3 log ■ 



27r 

as{s)CA 



271 



s 

7 



log^:^-^/?o^^°log- 



s 
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(26) 
(27) 



where = y|C'a — ^TpUj with Ca = 3 and Tp = \. These functions describe the 

total contribution for the emission of virtual particles (i.e. z = 1), with all invariants 
large compared to the cutoff /x, to the densities q{z,t) and g{z,t). The normalization 
is on the level of the cross section. For the invariant matrix element we thus find at 
the subleading level: 



M{pi,...,pn,gs,fj') = M{pi,...,pn,gs{s)) X 



1 1 9 

^ i=i ^ i=\ 



{2i 
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with Uq + Hg = n, and 



W^s,,^) = ^(log^A_31ogAj ,29, 

Vn..^) ^ H^(.og^^-i-,f".og-l) (30) 

Again we note that the running couphng notation in the Born-amphtude of Eq. (pHf ) 
denotes the renormahzation corrections of the Born amphtude and higher order cor- 
rections. The functions W^, correspond to the probabihty of emitting a virtual 
soft and/or coUinear gauge boson from the particle q, g subject to the infrared cutoff 
fi. Typical diagrams contributing to Eq. (|28| ) in a covariant gauge are depicted in Fig. 
^. In massless QCD there is no need for the label Wj or Wf, however, we write it for 
later convenience. The universality of the splitting functions is crucial in obtaining the 
above result. 



2.4 Renormalization group equation 



The solution presented in Eq. determines the evolution of the virtual scattering 
amplitude A1(pi, ...,p„, (7s, /i) for large energies at fixed angles and subject to the in- 
frared regulator fi. In the massless case there is a one to one correspondence between 
the high energy limit and the infrared limit as only the ratio s//i^ enters as a dimen- 



sionless variable fl^, |18|. Thus, we can generalize the Altarelli-Parisi equations (|24D 
and ( p5D to the invariant matrix element in the language of the renormalization group. 
For this purpose, we define the infrared singular (logarithmic) anomalous dimensions 



r,(t) . ; T,m . (31) 



Infrared divergent anomalous dimensions have been derived in the context of renor- 
malization properties of gauge invariant Wilson loop functionals |T9. In this context 



they are related to undifferentiable cusps of the path integration and the cusp angle 
PjPi/ gives rise to the logarithmic nature of the anomalous dimension. In case we use 



off-shell amplitudes, one also has contributions from end points of the integration |19 



The leading terms in the equation below have also been discussed in Refs. pO], |21 



and 1^ in the context of QCD. With these notations we find that Eq. ( |28D satisfies 

xM{pi,...,Pn,gs,i^) = (32) 

to the order we are working here and where Ai{pi, ...,pn, Qs, fJ') is taken on the mass- 
shell. The difference in the sign of the derivative term compared to Eq. (H) is due 
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Fi gure 3: The two counterterms contributing to the quark anomalous dimension — 
Qi^^j^-i {—Sqg + ^2)- Here JI denotes the MS dimensional regularization mass parameter. 
Due to divergences in loop corrections there are scaling violations also in the massless 
theory. 



to the fact that instead of differentiating with respect to log/i we use logs//i . The 



quark- antiquark operator anomalous dimension 7^5 = —Cp^^ enters even for massless 
theories as the quark antiquark operator leads to scaling violations through loop effects 
since the quark masslessness is not protected by gauge invariance and a dimensionful 
infrared cutoff needs to be introduced. Thus, although the Lagrangian contains no mqq 
term, quantum corrections lead to the anomalous scaling violations in the form of 7^^. 
The factor | occurs since we write Eq. ( P^ in terms of each external line separately^. 
For the gluon, the scaling violation due to the infrared cutoff are manifest in terms of an 
anomalous dimension proportional to the /5-function since the gluon mass is protected 
by gauge invariance from loop corrections. Thus, in the bosonic sector the subleading 
terms correspond effectively to a scale change of the coupling. Fig. || illustrates the 
corrections to the external quark- antiquark lines from loop effects. 

Except for the infrared singular anomalous dimension (Eq. all other terms 

in Eq. (521) are the standard contributions to the renormalization group equation for 



S-matrix elements In QCD, observables with infrared singular anomalous dimen- 
sions, regulated with a fictitious gluon mass, are ill defined due to the masslessness 
of gluons. In the electroweak theory, however, we can legitimately investigate only 
virtual corrections since the gauge bosons will require a mass. Eq. (^) will thus be 
very useful in the following sections. 



3 Logarithmic corrections in broken gauge theories 

In the following we will apply the results obtained in the previous sections to the case 
of spontaneously broken gauge theories. It will be necessary, at least at the subleading 

^In case of a massive theory, we could, for instance avoid the anomalous dimension term by 
adopting the pole mass definition. In this case, however, we would obtain terms in the wave function 
renormalization, and in any case, the one to one correspondence between UV and IR scaling, crucial 
for the vahdity of Eq. ([3^), is violated. 
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level, to distinguish between transverse and longitudinal degrees of freedom. The 
physical motivation in this approach is that for very large energies, s ^ = M^, 
the electroweak theory is in the unbroken phase, with an exact SU{2) x U{1) gauge 
symmetry. We will calculate the corrections to this theory and use the high energy 
solution as a matching condition for the regime for values of /x < M. 

We begin by considering some simple kinematic arguments for massive vector 
bosons. A vector boson at rest has momentum k'^ = (M, 0, 0, 0) and a polarization 
vector that is a linear combination of the three orthogonal unit vectors 

ei = (0,l,0,0) , ei = (0,0,l,0) , 63^(0,0,0,1) . (33) 

After boosting this particle along the 3-axis, its momentum will be k'^ = {Ek, 0, 0, k). 
The three possible polarization vectors are now still satisfying: 

k^e'^j =0 , e^j = -l . (34) 

Two of these vectors correspond to Ci and 62 and describe the transverse polarizations. 
The third vector satisfying (|3^ is the longitudinal polarization vector 

el{k) = {k/M,0,0,Ek/M) (35) 

i.e. e^(fc) = k'^ /M + 0{M/Ek) for large energies. These considerations illustrate that 
the transversely polarized degrees of freedom at high energies are related to the massless 
theory, while the longitudinal degrees of freedom need to be considered separately. 

Another manifestation of the different high energy nature of the two polarization 
states is contained in the Goldstone boson equivalence theorem. It states that the 
unphysical Goldstone boson that is "eaten up" by a massive gauge boson still controls 



its high energy asymptotics. A more precise formulation is given below in section |3T2 . 

Thus we can legitimately use the results obtained in the massless non-Abelian 
theory if we restrict ourselves to the transverse degrees of freedom at high energies. 
We will, however, show that to DL accuracy the results of Ref. DH can be used in 
connection with the Goldstone boson equivalence theorem. 

Another difference to the situation in an unbroken non-Abelian theory is the mixing 
of the physical fields with the fields in the unbroken phase. These complications are 
especially relevant for the Z-boson and the photon. 



3.1 Results for transverse degrees of freedom 

The results we obtain in this section are generally valid for spontaneously broken gauge 
theories, however, for definiteness we discuss only the electroweak Standard Model. 
The physical gauge bosons are thus a massless photon (described by the field A^) and 
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massive and Z bosons (described correspondingly by fields and Z,^).: 

= -l=(w^±tW^) (36) 



V2 

Z^ = cosO^Wj^ + sm0^B^ (37) 
= -sme„Wj^ + cos 9„B^ (38) 

Thus, amplitudes containing physical fields will correspond to a linear combination of 
the massless fields in the unbroken phase. The situation is illustrated schematically for 
a single gauge boson external leg in Fig. ^ In case of the bosons, the corrections 
factorize with respect to the physical amplitude. 

To logarithmic accuracy, all masses can be set equal: 

Mz ~ Mw ~ MHiggs ~ M 

and the energy considered to be much larger, ^/s ^ M. The physical fields are 
given in terms of the unbroken fields according to Eqs. (|5B|), ( P7| ) and (PB]). The 
left and right handed fermions are correspondingly doublets (T = 1/2) and singlets 
(T = 0) of the SU (2) weak isospin group and have hypercharge Y related to the electric 
charge Q, measured in units of the proton charge, by the Gell-Mann-Nishijima formula 
Q = T^ + Y/2. 

The value for the infrared cutoff fi can be chosen in two different regimes: 1) 
a/s ^ fi ^ M and 2) /i <C M. The second case is universal in the sense that it does 
not depend on the details of the electroweak theory and will be discussed below. In 
the first region we can neglect spontaneous symmetry breaking effects (in particular 
gauge boson masses) and consider the theory with fields B^, and W^^. One could of 
course also calculate everything in terms of the physical fields, however, we emphasize 
again that in this case we need to consider the photon also in region 1). The omission 
of the photon would lead to the violation of gauge invariance since the photon contains 
a mixture of the B^ and Wj^ fields. 

In region 1), the renormalization group equation (or generalized infrared evolution 
equation) ( ^2l) in the case of all rrii < M reads| 



x7W^(pi,...,p„,^,^',/x) = (39) 




^We exclude here top- Yukawa couplings which couple proportional to since they don't have an 
analogue in QCD. It is, however, not unlikely that those terms can also be included in the splitting 
functions fulfilling Altarelli-Parisi equations. Note also, that the amplitude on the right hand side is 
in general a linear combination of fields in the unbroken phase according to Eqs. (^6|), ( ^7| ) and (^Sf). 
In addition, in the electroweak theory matching will be required at the scale M and often on-shell 
renormalization of the couplings e and sin 0„ is used. In this case one has additional complications in 
the running coupling terms due to the different mass scales involved below M . Details are presented 
in section ^. 
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Figure 4: The schematic corrections to external gauge boson emissions in terms of the 
fields in the unbroken phase of the electroweak theory. There are no mixing terms 
between the Wj^ and B,, fields for massless fermions. We denote cos 9^ by c„ and 
sin6'w by Sw For final states, the corrections factorize with respect to the physical 
amplitude. In general, one has to sum over all fields of the unbroken theory with each 
amplitude being multiplied by the respective mixing coefficient. 
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where the index _L indicates that we consider only Ug transversely polarized external 
gauge bosons and n-w + '^b = ''^g- The two /3-functions are given by: 



P'ig'i-P')) 

with the one-loop terms given by: 

1 



dlogjf 
9 log 71^ 



-A 



-A 



' 87r2 



11 
— ( 
12 



Ao — TT^^A ~ g^gen 



24 



24^^^ 



(40) 
(41) 

(42) 



where n^en denotes the number of fermion generations ||2^, ^ and ri/j the number of 
Higgs doublets. The infrared singular anomalous dimensions read 



(43) 



where Tj and Yi are the total weak isospin and hypercharge respectively of the particle 
emitting the soft and coUinear gauge boson. Analogously, 



7, 



<?<? 



(44) 



The initial condition for Eq. (^) is given by the requirement that for the infrared 
cutoff /i^ = s we obtain the Born amplitude. The solution of (^) is thus given by 

M^ipU -^Pn, g, g', /i) = A^Born(Pl> --^Pn, g{s),g'{s)) 



X exp 



1 "9 



a{s) ^ ,^ , a' is) ( 
'Ti(T, + l) + 



47r 
ays' 



4tt 



log^ 



a is 



271 



Po + riB— — Ao log^ 



2tt 



1 "/ 



47r 



47r 



log ^ - 3 log — 



(45) 



where nw and ne denote the number of external W and -B fields respectively. The 
SU (2) X f/(l) group factors in the exponential can be written in terms of the parameters 
of the broken theory as follows: 



T.m + 1) + / ( y ) = + 9' (T.(r, + 1) - (T^r) + {Tf - sin' 6^ 



g' 



cos^ 6*,, 
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where the three terms on the r.h.s. correspond to the contributions of the soft photon 
(interacting with the electric charge Cj = Qig sin Oi^), the and the Z bosons, re- 
spectively. Although we may rewrite solution (^) in terms of the parameters of the 
broken theory in the form of a product of three exponents corresponding to the ex- 
changes of photons, and Z bosons, it would be wrong to identify the contributions 
of the diagrams without virtual photons with this expression for the particular case 
ef = 0. This becomes evident when we note that if we were to omit photon lines then 
the result would depend on the choice of gauge, and therefore be unphysical. Only for 
= 0, where the photon coincides with the B gauge boson, would the identification 
of the ef term with the contribution of the diagrams with photons be correct. 

We now need to discuss the solution in the general case. In region 1) we calculated 
the scattering amplitude for the theory in the unbroken phase in the massless limit. 
Choosing the cutoff /i in region 2), /x -C M, we have to only consider the photon 
contribution. In this region we cannot necessarily neglect all mass terms, so we need 
to discuss the subleading terms for QED with mass effects. If rrii -C /i, the results from 
the massless QCD discussion of section |2.3| can be used directly by using the Abelian 
limit Cp = 1. In case fi <^ rrii we must use the well known next to leading order QED 
results, e.g. |26], and the virtual probabilities take the following form for fermions: 



' (^(W#-31og^ 

log^-l)21ogf 

+ W^-31og^ 



(47r)2 



ITLi ^ H 



(46) 



Note, that in the last equation the full subleading collinear logarithmic term is 
used in distinction to Ref. 



26 



In the explicit two loop calculation presented in Ref. 
it can be seen that the full collinear term also exponentiates at the subleading 
level in massive QED. For bosons we have analogously: 



(47r)^ 



(47) 



In addition we have collinear terms for external on-shell photon lines f\ from fermions 
with mass rrij and electromagnetic charge Cj up to scale M: 



3 47r2^^C^Og ^2 
1 y-"/ A. Mi lo(T Ml 



/i <^ rrtj 



(4J 



Note that automatically, w]{M^, M^) = 0. At one loop order, this contribution cancels 
against terms from the renormalization of the QED coupling up to scale M. For 

thank the authors of Ref. for clarifying this point. 
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external Z-bosons, however, there are no such colhnear terms since the mass is large 
compared to the m^. Thus, the corresponding RG- logarithms up to scale M remain 
uncanceled. 

The appropriate initial condition is given by Eq. (|45| ) evaluated at the matching 
point /i = M. Thus we find for the general solution in region 2): 



M^iPi,...,Pn,g,g',f^) = Mi^^^{pi,...,Pn,gis),g'{s)) 



i=l 



a(s'' 



Ti{Ti + l) + 



a [s J 
47r 



log' 



M2 



a is) 



+ nw^z—Po + riB— — P'o log — 



27r 



1 "/ 

^ k=l 



2n 
11 



47r 



47r 



2 



X exp 



1 1 



i=\ 



= M 

X exp 



-^E«^7(M^/i') 



Born(Pl.--.Pn,^(s),^'(s)) 



1 "9 

2 ^ 



47r 



27r 



/5o + 



a [s 



An 
-P'o] lo^ 



2 



S 

M2 



1 "/ 

2 ^ 



-i^un + 1) + 



X exp 



1 " 

-2g' 



,(4vr) 



2 log 



47r 



log' -7-pj — 3 log — — r 



M2 



M2 



-log— + 21og — log^ + 31og— I 

rriiM mf mf fi'^ 



mfU l"Z,le? , M2 



(49) 



The last equality holds for /i ^ mj ^ M and we have absorbed all /5-function terms 
not related to external lines into redefinitions of the scales of the couplings. It is 
important to note again that, unlike the situation in QCD, in the electroweak theory 
we have in general different mass scales determining the running of the couplings of the 
physical on-shell renormalization scheme quantities. We have written the above result 
in such a way that it holds for arbitrary chiral fermions and transversely polarized 
gauge bosons. In order to include physical external photon states in the on-shell 
scheme, the renormalization condition is given by the requirement that the physical 
photon does not mix with the Z-boson. This leads to the condition that the Weinberg 
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rotations in Fig. ^ at one loop receive no RG-corrections. Thus, above the scale M 
the subleading collinear and RG-corrections cancel for physical photon and Z-boson 
states. For physical observables, soft real photon emission must be taken into account 
in an inclusive (or semi inclusive) way and the parameter /i^ in (^9]) will be replaced by 
parameters depending on the experimental requirements. This will be briefly discussed 
in the following section. 



3.1.1 Semi inclusive cross sections 

In order to make predictions for observable cross sections, the unphysical infrared cutoff 
jj?' has to be replaced with a cutoff /ig^p, related to the lower bound of fc^ of the other 
virtual particles of those gauge bosons emitted in the process which are not included 
in the cross section. We assume that /ig^.^ < M^, so that the non-Abelian component 
of the photon is not essential. The case ^1^^ > is much more complicated and is 
discussed in Ref. [0 through two loops at the DL level. 

We again only discuss transversely polarized external gauge boson in the Born 
process and can write the expression for the semi-inclusive cross section: 



da^ipi,. . .,Pn,g,g',fiexp) = da^^^^^-^{pi, . . . ,pn, g{s), g'{s), fi) 

X exp(wj^p(s, rui, /i, /ie^.p)) 



(50) 



In the soft photon approximation we have: 



e. 

o2 



(47r)2 ^j=l 

-2 log 



■log' 

log;; 



+ log' 



Slog 



1 21og^ + log 



2 s 



log A - 1 



(51) 



where the upper case applies only to fermions since for we have fi < M in region 
2). Since the upper bound on of the photons which are allowed to be radiated is 
less than M^, we must use the cut-off /i^ < and, consequently, (^) for the matrix 
element of the non-radiative process. Therefore, we obtain* 

da^ipi, ...,Pn,g, g', f^exp) = C?C^Born(Pl> • • • ,Pn, g{s),g'{s)) 



X exp <^ - 2 



i=l 



^ Air 

ais" 



a {sj 

4:71 



1 2 S 



a'(s) \ s 
nw^Po + ns-^f3'o 1 log ^ 



®The notation here is again simplified in the sense that for Z-boson and 7 final states one has to 
include the mixing correctly as described above. 
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X exp 



X exp 



i=l 



i=l 



2 2\ 



i=l 
^2 



E 



i=l 



1 2 -5 

+ log — 

mf 



(52) 



where we use fi -C mj. The dependence in this expression cancels and the semi- 
inclusive cross section depends only on the parameters of the experimental require- 
ments. 

Eq. ( |52D contains all leading double and single logarithms to cross sections^ con- 
taining arbitrary numbers of external fermions and transversely polarized gauge bosons. 
We have only assumed that all masses are not larger than the electroweak scale M and 
impose a cut on the the allowed values of emitted real gauge bosons < nl^p < M^, 
i.e. up to the weak scale we only need to consider real QED effects. 



3.2 Longitudinal degrees of freedom 

In this section we discuss if results obtained from the massless unbroken phase of the 
SU{2) X U{1) theory, where due to gauge invariance we have only transverse physical 
degrees of freedom, can be extended to the full theory including longitudinal vector 
bosons. This point of discussion is necessary and important since the longitudinal 
degrees of freedom don't decouple at high energies and could give crucial clues to 



potentially strong dynamical effects for large Higgs masses mn ~ ITeV ^4 . 

The connection between the strategy pursued for the transverse degrees of free- 
dom and the corrections to longitudinally polarized vector bosons at high energies is 
provided by the Goldstone boson equivalence theorem |2^. It states that at tree level 
for S-matrix elements for longitudinal bosons at the high energy limit M^/ s — > can 
be expressed through matrix elements involving their associated would be Goldstone 
bosons. We write schematically in case of a single gauge boson: 

A<(W^f,V'phys) = A^(0^,Vys)+(^(^) (53) 



^We emphasize again that we did not consider angular logarithms which can be sizable and should 
be calculated at least to one loop order. 
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M{Zl, V'phys 



-M(0,^phys) + O 



(54) 



The problem with this statement of the equivalence theorem is that it holds only at tree 
level p9| , 0. For calculations at higher orders, additional terms enter which change 
Eqs. (H) and (H). 

Because of the gauge invariance of the physical theory and the associated BRST 
invariance, a modified version of Eqs. 



rA^(Ty,^(fc),Vys) 



31) and (|5J) can be derived p9| which reads 
± 



C^M^M{(j)^{k),iljp^y,) + 
C,M,M{(p{k),ijpi,ys 



(55) 
(56) 



where the multiplicative factors and depend only on wave function renormal- 
ization constants and mass counterterms. Thus, using the form of the longitudinal 
polarization vector of Eq. (pSj) we can write 



7U(W±(A;),Vys) 
M{ZL{k),iljphys) 



C„A<(0±(A;),^phy.) + O 
C,A^(0(A;),Vys) + (^( 



AC 



(57) 
(58) 



Thus we see that in principle, there are logarithmic loop corrections to the tree level 
equivalence theorem^. In addition, for longitudinal gauge bosons we also have log- 
arithmic corrections with Yukawa terms [Q. On the one hand, this means that the 
method of section ^A] should be used with caution to obtain the relevant subleading 
terms. Thus we should consider these corrections separately. On the other, since 
the corrections are at most logarithmic, it means that the results of Ref. can be 
extended to the longitudinal sector as well. Thus we fincQ for fi^^p < M: 

da^Pi, . . .,Pn,g,g',fXexp) = da^^^^{pi, . . .,pn,g,g') 



E(^T.(T. + l) + ^(i^) ilog^-l, 



a 



X exp 
X exp 



1=1 

n 



-E(- 



47r ^ ' 47r 
s,/i2) -w^ 



M2 



DL,„ ,.2^ ..,DL,^2^ 



+ w 



DL 

exp 



i=l 



(59) 



An exception is the background field gauge where the Ward-identities guarantee that the factors 
Cw — 1 and Cz = 1 to all orders ||3l| . It should thus be investigated if subleading corrections can also 
be obtained from the Goldstone boson equivalence theorem. 

^^For longitudinally polarized Z-boson final states there are no mixing terms since the photon has 
only transverse polarization states. Thus one needs to only include the associated Goldstone boson (j) 
at the DL level. 
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Figure 5: The pictorial Goldstone boson equivalence theorem for ly-pair production 
in e^e~ collisions. The correct DL-asymptotics for longitudinally polarized bosons are 
obtained by using the quantum numbers of the charged would be Goldstone scalars at 
high energies. 

where the index || indicates the cross section for longitudinally polarized gauge bosons, 
while the field indicates that the appropriate fields and quantum numbers on the 
r.h.s. in Eq. (|59|) are those of the associated would be Goldstone bosons. 

Thus, we have shown that all DL corrections can be summed to all orders by 
employing the evolution equation approach of Ref. |l[ in connection with the Goldstone 
boson equivalence theorem. 



4 Comparison with explicit results 

In this section we compare our results obtained in the previous sections with known 
results in special cases and one loop calculations. In Ref. QCD-results for the 

Sudakov form factor were generalized to the high energy electroweak theory^^. Since 
the general strategy pursued is the same as in Ref. we of course agree with their 
result for leading and subleading electroweak corrections to e^e~ — > ff to all orders. 

A very important check is provided by the explicit one-loop corrections of Ref. 
for high energy on-shell W^-pair production in the soft photon approximation. In the 
following, the lower index on the cross section indicates the helicity of the electron, 
where el denotes the left handed electron. We summarize the relevant results for 
e^el — > W^Wl[, e^eZ — > W^^W^^ and e^e^ — > W^^W^~ for convenience as follows: 



da\ f da 



\ Born 



- 



dni _ 




l + 2c2 +8ct, , s , l-2ci + Act 



""^^In addition, angular terms at the one loop level were calculated which we do not consider in this 
work. 
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+3 log — + 2 log 



-.EQlAT^log-l 




(60) 



l-2cl + Act 



2cls 



2 



log^ 



M2 



s s 

log — + log 77^ 

mi 



10c: 



5 S 

log — + log 




(61) 



(62) 



where the last line in Eq. (|60| ) corresponds to a sum over all fermions contributing 
to the coupling renormalization (with multiplicity Nc = 3 for quarks and Nc = 1 
for leptons). This contribution can be included in the scale of the running on-shell 
charge aesiM"^) @]- For the longitudinal cross sections we are only concerned with 
DL corrections. The Born cross sections are given by: 



Born 




1 u' + f 



\_ 

647r2s 164ct 

1 



sin^e 



647r2s 4cl 



sin'^ 



(63) 
(64) 
(65) 



where we keep the angular dependence. In Eq. (^) a sum over the two transverse 

polarizations of the (++ and ) is implicit. These expressions demonstrate that 

the longitudinal cross sections in Eqs. (|64D and (|65|) are not suppressed with respect 



to Eq. (|63|). On the other hand, ^ is mass suppressed 0. 

Eqs. (0), (|6TD and (|62D were of course calculated in terms of the physical fields of 
the broken theory and in the on-shell scheme. We denote Cw = cos6'w and = sin6'„ 
respectively. In order to compare with the results of section ^ we listed the relevant 
quantum numbers in Table |I|. For comparison with Eq. (^) and to logarithmic 
accuracy, we can absorb the running coupling effects from our massless scheme to the 
on-shell scheme as follows. The Born cross section in our approach is proportional to 



22 



g'^ (see Eq. (|63|)). The coupling renormalization above the scale M is given by: 



i-(M^) f 1 _ im flic, - 1„, 



n 



gen 



Below the scale where non-Abelian effects enter, the running is only due to the elec- 

_„2 



tromagnetic coupling and we write g'^{M'^ 



with 



1 p2 • 



(67) 



We therefore observe that the running coupling terms proportional to log cancel 
for this process with the subleading contributions from the virtual splitting functions 
(see Eq. (|49|) ) and what remains are just the Abelian terms up to scale M. Thus for 
Eq. (1^01) we obtain from Eq. (K^) at the one loop level: 
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H 

87r2 



log 



2 e2 ^^2,,,, M2 
+ 3^EQ|iV^log-, 



(68) 



Eq. (|68D agrees with Eq. (|60D, which are both valid in the soft photon approximation. 
Here and below we assume that AE < M and fiexp < M. Analogously in the DL 
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T 


Y 





p 


1/2 











-2 


-1 




1/2 


1 


1 


ei 





2 


1 




1 





±1 


0± 


1/2 


±1 


±1 



Table 1: The quantum numbers of various particles in the electroweak theory. The 
indices indicate the helicity of the electrons. We neglect all mass terms. For longitudi- 
nally polarized gauge bosons, the associated scalar Goldstone bosons describe the DL 
asymptotics. 

approximation, it is straightforward to check the validity of our results for Eqs. (|6TD 
and emphasizing again that in this case we need to use the quantum numbers of 
the associated Goldstone bosons, see Fig. ^ and Tab. 0. 

Thus we have verified that our results, calculated in terms of the unbroken massless 
fields, give the correct leading and subleading logarithms in transversely polarized Im- 
pair production at the one loop level. For longitudinally polarized VT-pairs the correct 
DL asymptotics is reproduced. 

5 Conclusions 

In this paper we considered the calculation of virtual next to leading electroweak 
corrections at energies much larger than the electroweak scale when all particle masses 
can be neglected. We follow the same approach as in Ref. which consists of using 
the fields of the unbroken theory to obtain logarithmic corrections with the infrared 
evolution equation method in different regions of the infrared cutoff. When particle 
masses can be neglected there is a one to one correspondence between the high and 
low energy scaling behavior and the evolution equation can be formulated in terms of 
the renormalization group with infrared singular anomalous dimension. The next to 
leading kernel can then be obtained from the virtual contribution to the Altarelli-Parisi 
splitting functions. 

For external gauge boson emission one can use the above approach for transverse 
degrees of freedom. For fermions and external states, the next to leading correc- 



24 



tions exponentiate with respect to the physical Born amphtudes. For Z boson and 7 
final states, one needs to include the effect of mixing appropriately. For these final 
states we have exponentiation with respect to the amplitudes containing the fields of 
the unbroken theory but not with respect to the physical Born amplitude. 

For longitudinal degrees of freedom, one can use the Goldstone-boson equivalence 
theorem to obtain the correct DL asymptotic behavior. These terms are found to ex- 
ponentiate as well. Loop corrections, however, lead to additional corrections including 
Yukawa terms at the subleading level and should be considered separately. 

We also compared our results with exact one loop calculations in the high energy 
approximation and could reproduce the leading and subleading terms for transversely 
polarized VF-pair production in e^e~ collisions and the DL corrections for the longitu- 
dinal degrees of freedom. 

Finally we note that there are of course terms which we have neglected in this 
analysis. As mentioned above, there are angular logarithms of the form log -j^la log ^, 
which in general could be significant and should be computed separately. We also 
omitted all mass-logarithms of the form log log and top- Yukawa terms. For 
the latter, it might be possible to include them consistently into the virtual splitting 
functions. For longitudinal degrees of freedom it would also be very helpful to calculate 
subleading corrections. 

In conclusion, for future high energy experiments in the multi-TeV energy regime, 
the leading high energy behavior of general scattering amplitudes can be an important 
ingredient to study the effect of new physics expected in precisely this range. 
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